On the applicability of bosonization and the Anderson- Yuval methods at the 
strong-coupling limit of quantum impurity problems 
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The applicability of bosonization and the Anderson- Yuval (AY) approach at strong coupling is 
investigated by considering two generic impurity models: the interacting resonant-level model and 
the anisotropic Kondo model. The two methods differ in the renormalization of the conduction- 
electron density of states (DoS) near the impurity site. Reduction of the DoS, absent in bosonization 
but accounted for in the AY approach, is shown to be vital in some models yet superfluous in others. 
The criterion being the stability of the strong-coupling fixed point. Renormalization of the DoS is 
essential for an unstable fixed point, but superfluous when a decoupled entity with local dynamics 
is formed. This rule can be used to boost the accuracy of both methods at strong coupling. 

PACS numbers: 72.10.Fk, 72.15.Qm, 73.63.Kv 
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— Introduction. Several classic models in condensed- 
matter physics show logarithmic behavior at high ener- 
gies, followed by qualitatively different behavior at low 
energies. Notable examples include the x-ray absorp- 
tion problem^ the Kondo Hamiltonian^ the interacting 
resonant-level model (IRLM) j&i and different variants of 
two- level systems (TLS). 5,6 Historically devised to model 
real impurities in bulk samples, many of these Hamiltoni- 
ans have recently found new realizations and generaliza- 
tions in quantum dots and other confined nanostructures. 

A distinguished place in the theory of such quantum 
impurities is reserved to Abelian bosonization 7 " and the 
Anderson- Yuval (AY) approach;^ which remain among 
the most powerful and versatile analytical tools in this 
realm. With numerous applications over the last forty 
years, it is surprising that the applicability of neither ap- 
proach has ever been studied systematically for strong 
couplings. In bosonization, the bare couplings are gener- 
ally assumed to be weak. Strong static interactions are 
often included ad hoc in terms of their scattering phase 
shift. The AY method, which maps the original impurity 
problem onto an effective Coulomb gas, is presumably 
nonperturbative in certain couplings. However, it typi- 
cally fails to reproduce the correct scaling equations even 
at the next-to-leading order. 4,10 A reliable extension of 
these approaches to strong couplings is highly desirable. 

The goal of the present paper is to critically test the 
accuracy of these leading analytic methods away from 
weak coupling, and to propose an operational extension 
to strong couplings. To this end, we resort to Wilson's 
numerical renormalization group 11 (NRG), and to two 
generic classes of models as test beds: the IRLM and the 
anisotropic Kondo model. Our analysis highlights the 
role of the reduction in the conduction-electron density 
of states (DoS) near the impurity site, which may hinder 
the efficiency of other essential couplings (e.g., tunnel- 
ing in the IRLM). This reduction of the DoS, absent in 
bosonization but included in the AY approach, proves 



vital in some models and superfluous in others. It is es- 
sential in cases where the strong-coupling fixed point is 
unstable, but superfluous in models where a decoupled 
entity with local dynamics is formed at strong coupling. 
Hence, the accuracy of bosonization and the AY approach 
can be significantly enhanced by selectively incorporat- 
ing the DoS renormalization factor to match the case in 
question. 

The reduction of the local conduction-electron DoS is 
best seen for a simple model where electrons scatter elas- 
tically off a point-like impurity (s-wave scattering). The 
renormalized DoS at the impurity site takes the form 12 



q(lu « Ep) = go cos 2 8 



(1) 



where go is the unperturbed DoS, Ep is the Fermi en- 
ergy, and 5 is the scattering phase shift. Since S — > tt/2 
for resonant scattering, this implies g(u> ~ Ef) — > 0. 
This fact may have a dramatic effect, as exemplified be- 
low by the two-channel IRLM. A strong local Coulomb 
repulsion suppresses the DoS at the vicinity of the impu- 
rity, reducing the hopping rate between the impurity and 
the bands. Since reduction of the DoS is independent of 
the interaction sign, it equally applies to an alternating 
potential. The case of a TLS with a single coupling (the 
commutative model 5,6 ) is qualitatively similar. 

One may expect the same to occur in the anisotropic 
Kondo model or the non-commutative TLS with electron- 
assisted hopping. For example, consider the single- 
channel Kondo model (1CKM) with a large XXZ 
anisotropy: J z ^> \J±\, with J x — J y = J±. In the 
spirit of the AY philosophy^ one may first treat the larger 
coupling J 2 before incorporating the smaller J±. In the 
absence of Jx, a large J z reduces the local DoS at the 
impurity site independent of the orientation of the im- 
purity spin. Incorporating J± at the next step, its effi- 
ciency is expected to be hindered by the reduced DoS, to 
the extent that it diminishes in the limit J z — > oo [when 
5 — > 7r/2 and g(u> « Ep) — > 0]. Surprisingly, this is not 
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what we find with the NRG. Rather, spin flips remain 
governed at large J z by the bare transverse coupling Jj_ . 

To unravel the governing rule, we conduct a detailed 
comparison between Wilson's NRG, bosonization, and 
the AY method, applied separately to the multichannel 
Kondo and IRLM models. Applicability of the latter two 
approaches at strong coupling is shown to depend cru- 
cially on the stability of the strong-coupling limit. When- 
ever a decoupled entity with local dynamics is formed 
(i.e., a stable strong-coupling fixed point is reached), 
then the DoS renormalization factor is superfluous and 
bosonization works well. If, however, the strong-coupling 
limit is unstable, then the DoS renormalization factor is 
essential and the AY approach works well. The above 
classification pertains to non-commutative models. For 
commutative couplings the AY method always applies as 
one can always reorder the perturbation series. 

Prompted by these findings we proceed to re-examine 
the "intimate relation" between the IRLM and the 
anisotropic lCKMi 4 - Close correspondence is established 
between the models in case of the single-channel IRLM, 
but not in the case of multiple screening channels. 

— Interacting resonant-level model. In the IRLM^i^ a 
ID electron gas is coupled to a spinless impurity level 
by two distinct mechanisms: a hopping matrix element 
V and a short-range Coulomb repulsion U . The hop- 
ping rate is enhanced for weak repulsion, but is gen- 
erally suppressed at large U due to a reduction in the 
conduction-electron overlap integrals between a vacant 
level and an occupied on o 13 > 14 (the so-called orthogonal- 
ity catastrophe 1 ^) . Consequently, the hopping rate tends 
to develop a maximum at some intermediate coupling U, 
whose value is pushed toward weak coupling as the num- 
ber of screening bands N is increased*^ This behavior 
stems from an enhancement of the orthogonality effect 
with increasing N. 

Interest in the IRLM has been recently rekindled by a 
Bethe Ansatz solution of a two-lead version of the model 
under nonequilibrium conditions. 16 In its multichannel 
form, the Hamiltonian reads H = Ho + Hi + TL2, with 

jV-l 

Ho = v F{k-k F )a{ n a kn + t d d' i d , (2) 

n=0 0<fc<2fc F 

Hi = Uj^ (4a„-^)(^-£) , (3) 
H 2 = V(d)a + ald} . (4) 

Here, a\ n creates an electron with momentum k in the 
nth band, d) creates an electron on the level, kp and 
Vf are the Fermi momentum and Fermi velocity, respec- 
tively, €d is the level energy, U is the Coulomb repulsion, 
and V is the tunneling amplitude into the n = band. 
The operator a) n = (1/ VJ7) J^k a ln > where Af is the num- 
ber of lattice sites, creates a localized band electron at 
the impurity site. Note that H is particle-hole symmetric 
for Ed = 0, the case of interest here. 
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FIG. 1: (Color online) Renormalized level width for the IRLM 
with up to 4 screening channels, as obtained by the NRG, 
bosonization, and the AY approach. Here To = ttqoV 2 with 
V/D = 0.02 (we use g D = 1/2). The AY approach (panel 
a) works quite well for N = 2,3,4, but fails for N = 1. 
Bosonization (panel b) works well for N = 1,3,4, but pre- 
dicts an exact mapping^ between U — *■ and U — > 00 for 
N — 2, and thus a saturated width. Note that the AY ap- 
proach systematically underestimates r rcn at large U whereas 
the opposite is true of bosonization. 

We study the IRLM using Wilson's NRG, bosonization 
and the AY approach. Since bosonization and the NRG 
are frequently used, we refer the reader to Refs. and [Til 
for details of these methods. In the following we briefly 
review the AY approach, which relies on a mapping of 
the impurity problem onto an effective ID Coulomb gas 
of multicomponent charges. The AY mapping is nonper- 
turbative in the Coulomb repulsion U, which determines 
the different charge components through its associated 
phase shift S = — arctan(7rpotV2). Here go is the bare 
conduction-electron DoS. The hopping amplitude V fixes 
the fugacity of the gas, which is given in turn by 

y = V(g To) 1/2 cosS. (5) 

Here To = 1/Dq is a short-time cutoff, with Do the bare 
bandwidth. The cos 8 that appears in Eq. © encodes the 
DoS renormalization. A similar mapping, only without 
the cos S, can be derived using Abelian bosonization. In- 
corporating U by means of its associated phase shifty 
an identical ID gas is obtained with y — V(qoTo) 1 / 2 . 

The Coulomb gas is next treated by progressively in- 
creasing the short-time cutoff while simultaneously renor- 
malizing the gas parameters so as to leave the partition 
function invariant. This results in renormalization-group 
(RG) equations for the parameters of the Coulomb gas 14 
which are perturbative in the fugacity y (namely, V) but 
nonperturbative in U. To illustrate the basic iterative 
step, suppose that the short-time cutoff has already been 
increased from its bare value To = I/-D0 to r > To. Fur- 
ther increasing the cutoff to r + St requires two opera- 
tions: (i) integration over charge pairs whose separation 
falls in the interval (r, r + St) , and (ii) rescaling of r by 
t + St. Consecutive charges, having opposite signs, leave 
no net charge behind. However, they do possess a dipole 
moment that acts to screen the interaction between the 
charges that remain. Integration over the close-by charge 
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pairs can therefore be absorbed into a renormalization of 
the remaining charges. On the other hand, the rescaling 
of t is absorbed into a renormalization of the fugacity y, 
as described by the following set of RG equations^ 



dy 
dlnr 

dlnr 



N-l 



2 (<5„o + z n ) y 2 



(6) 
(7) 



Here <5„o is the Kronecker delta, while the charge compo- 
nents z n take the bare value z = 25 /n. Contrary to usual 
dynamical scaling equations, the DoS is also modified in 
this procedure due to the rescaling of t. However, this 
difference is only formal. Either strategy can be pursued. 

Equation ([6]) pertains to the fugacity y. It can equally 
be written as a scaling equation for the level width 
r = iry 2 /r, which serves as the low-energy cutoff in 
the problem. Specifically, the perturbative RG procedure 
terminates at 1/r ~ T, when the fugacity y becomes of 
order 1. Whether this condition is met or not depends 
on the values of N and S. To see this, consider a suffi- 
ciently small i/o such that the renormalizations of z n can 
be ignored. Equation ([6]) then becomes 



= -(l-2z- Nz 2 ) y. 
dlnr 2 V > y 



(8) 



Whether y is relevant or not depends on the sign of the 
expression in the brackets. Since — 1 < z < for repulsive 
interactions, y is always relevant for N < 3. However, it 
turns irrelevant for N > 3 if U is made sufficiently large. 
The system flows then to a decoupled level. Careful anal- 
ysis of the transition between a strongly coupled and a 
decoupled level shows that it is of the Kosterlitz-Thouless 
type?i£ analogous to the ferromagnetic-antiferromagnetic 
transition line of the anisotropic Kondo model. Impor- 
tantly, bosonization and the AY approach predict the 
same critical coupling U c as V — > 0. 

Solution of Eq.© in the regime where y is relevant 
yields the renormalized level width, or cutoff scale, 



r re 



D o y 



2/(l-2z-Nz 2 ) 



(9) 



Here yo is the bare fugacity of Eq. ([5]) . For either N = 1 
or N = 2, one can substitute z ~ — 1 in Eq.© to obtain 

2/(3 — N) 

r rcn ~ Doy at large goU. Hence r ron is strongly 

suppressed as go U — > oo due to the cos<5 that appears 
in yo- By contrast, r ren saturates in bosonization, where 
the DoS renormalization factor is absent. 

Figure [T] compares the renormalized level width r ron of 
the multichannel IRLM, as obtained by our three meth- 
ods of interest. Within the NRG, r rcn was defined from 
the T — ► charge susceptibility of the level according 
to r ron = l/7r% c . In the AY approach and bosoniza- 
tion, r ren was obtained from a full numerical solution of 
Eqs.© and 0, with and without the cos<5 in Eq.©. 

While both the AY method and bosonization work 
quite well for N > 2, only the former approach succeeds 
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FIG. 2: (Color online) The Kondo temperature of the one- 
and two-channel Kondo model as a function of J z for fixed 
QoJ± ~ 0.1. While bosonization works quite well for the 
1CKM, the AY approach incorrectly predicts a vanishing Tk 
as J z — » oo. The roles are reversed for the 2CKM. Here 
the AY method is qualitatively correct, whereas bosonization 
predicts^ an exact mapping between J z — > and J z — » oo, 
and thus a saturated Tk- The dotted-dashed line shows a 
one-parameter fit [the prefactor of Ax in Eq.(47) of Ref. |2l|] 
to the strong-coupling expansion of the 2CKM. 



in tracing the NRG for N = 2. Bosonization fails to pro- 
duce the suppression in r rcn at large U, which stems from 
the renormalized DoS. By contrast, the AY approach fails 
to generate the saturation in r ren for N — 1 and large 
U, which bosonization captures quite well. Hence, the 
DoS renormalization factor is superfluous in this case. 
The source of distinction between N = 1 and N = 2 
is nicely elucidated by a strong-coupling expansion^ in 
1/U. Whereas a decoupled entity with local dynamics is 
formed when N = 1, for N = 2 the strong-coupling fixed 
point is unstable. A renormalized IRLM is recovered^ 
with dynamics that depends on the renormalized DoS. 
Relevance of the DoS renormalization depends then on 
the stability of the strong-coupling fixed point. As shown 
below, the same criterion applies to the Kondo model. 

— Anisotropic single- channel Kondo model. The 
anisotropic 1CKM has been intensely studied over the 
year o 8 ' 9 ' 10 ! 19 as a paradigmatic example for strong corre- 
lations. It describes the spin-exchange interaction of an 
impurity spin S with the local conduction-electron spin- 
density s, as modeled by the Hamiltonian term 



H 



int 



(10) 



In the antiferromagnetic regime, J z > — |J_i_|, the system 
flows to the strong-coupling fixed point of the isotropic 
model regardless how large the anisotropy is. 

Similar to the hopping V in the IRLM, the transverse 
Kondo coupling J± is attached a factor of cos 2 5 with 
S = — arctan(7rpo</z/4) upon mapping the 1CKM onto an 
effective ID Coulomb gas using the AY approach. This 
factor, which stems from the form of the electronic Green 
function,^ is absent in bosonization, and is omitted in 
the original works of Anderson and collaborators^ Its 
inclusion has profound implications, as the effect of spin 
flips (and consequently the Kondo temperature) vanishes 
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FIG. 3: (Color online) Renormalized level width of the one- 
and two-channel IRLM vs the Kondo temperature of the 
1CKM, obtained using the NRG. Here J± = V&V. As 
predicted by bosonization, there is close correspondence be- 
tween the 1CKM and the N = 1 IRLM upon equating 
V2S Z + 7r(V2- l)/2 with «5u- 

in the limit 5 — ► \ (i-e., J z — > oo). If these considerations 
are correct, then the NRG should give the same result as 
J z — > oo, which turns out not to be the case. 

Figure fj] compares the Kondo temperature Tk ob- 
tained by our three methods of interest. Within the 
NRG, Tk was defined from the T — ► impurity spin sus- 
ceptibility according to Tk — l/^Xs- In the AY approach 
and bosonization, it followed from a full numerical solu- 
tion of the RG equations 9 , with and without the cos 2 5 fac- 
tor attached to J±. Evidently, bosonization works quite 
well for the 1CKM, reproducing the saturation of the 
Kondo temperature as J z — > oo. The AY prediction of 
a vanishing Tk is clearly discredited by the NRG, prov- 
ing the redundancy of the DoS renormalization factor. 
As anticipated, a decoupled entity is formed at large J z , 
signaling the stability of the strong-coupling fixed point. 

A critical test of our picture is provided by the 
anisotropic two-channel Kondo model (2CKM), whose 
strong-coupling fixed point is known to be unstable. 
Instead, the model flows to an intermediate-coupling, 
non- Fermi- liquid fixed point, characterized by anomalous 
thermodynamic and dynamic properties. 20 Similar to the 
two-channel IRLM, we expect the DoS renormalization 
to be essential in this case. The results shown in Fig[2] 
well support our picture. While bosonization predicts 2 ^ 
an exact mapping between J z — > and J z — > oo, and 
thus a saturated Tk, the AY approach correctly repro- 



duces the vanishing of Tk- Though quantitatively less 
accurate at intermediate J z , agreement with the NRG is 
clearly very good both at small and large coupling. 

Above two screening channels, the anisotropic Kondo 
model undergoes a Kosterlitz-Thouless transition with 
increasing J z > to a ferromagnetic-like stated Since 
spin flips are suppressed to zero, the distinction between 
bosonization and the AY approach looses its significance 
at strong coupling, similar to the IRLM with N > 3. 

— Comparison of the two models. Prompted by these 
results, we have set out to carefully test the accepted 
mapping^ of the one-channel IRLM onto the 1CKM, as 
the mapping involves large couplings. Within bosoniza- 
tion, one finds the following correspondence of parame- 
ters^ V <-> Jj_/VB and <5tj <-> V2S z +w(V2- l)/2, with 
<5tj = — arctan(7T0o?7/2) and S z = — arctan(7T£>o</z/4). 
Our NRG results for the low-energy scales of both mod- 
els are summarized in Fig|3J Evidently, there is close 
correspondence between the two models using the above 
mapping of parameters, confirming the predictions of 
bosonization. Note that Tk varies by a factor of 30 in 
FigEl The agreement does not extend to the two-channel 
IRLM, which similarly flows to a strong-coupling Fermi- 
liquid fixed point (unlike the non- Fermi-liquid fixed point 
of the 2CKM). The DoS renormalization factor, absent 
in the 1CKM, proves essential in this case. 

— Conclusions. We have critically examined the ac- 
curacy of the AY and bosonization methods away from 
weak coupling by considering two generic impurity mod- 
els. Reduction of the conduction-electron DoS, accounted 
for by the AY approach but absent in bosonization, was 
shown to be vital in the case of an unstable strong- 
coupling fixed point, yet superfluous in models where 
a decoupled entity with local dynamics is formed. The 
two methods thus display complementary accuracies at 
strong coupling, controlled by the stability of the strong- 
coupling fixed point. Accuracy of these powerful methods 
can thus be significantly enhanced by selectively incor- 
porating the DoS renormalization factor, making them 
adequate tools for tackling strong-coupling physics. 
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